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I. INTRODUCTION 



Understanding the origin of the patterns of neutrino masses and mixing, emerging from 
the neutrino oscillation, 3 H (3— decay, etc. data is one of the most challenging problems in 
neutrino physics. It is part of the more general fundamental problem in particle physics of 
understanding the origins of flavour, i.e., of the patterns of the quark, charged lepton and 
neutrino masses and of the quark and lepton mixing. 

At present we have compelling evidence for existence of mixing of three light massive 
neutrinos z/j, % = 1,2,3, in the weak charged lepton current (see, e.g., jl|). The masses rrii 
of the three light neutrinos z/j do not exceed approximately 1 eV, m; ^ 1 eV, i.e., they are 
much smaller than the masses of the charged leptons and quarks. The three light neutrino 
mixing is described (to a good approximation) by the Pontecorvo, Maki, Nakagawa, Sakata 
PMNS) 3x3 unitary mixing matrix, £/pmns- 111 the widely used standard parametrisation 
l|j t^PMNS is expressed in terms of the solar, atmospheric and reactor neutrino mixing angles 
#12) #23 and # 13 , respectively, and one Dirac - 5, and two Majorana [2] - a 2 \ and a 31 , CP 
violating phases: 

Upmns = U = 1/(6*12, 6» 23 , 6>i3, 5) Q{a 21 , a 3 i) , (1) 



where 
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and we have used the standard notation = cos^-, = sin^y, and 



(2) 



Q = diag(l, e 



i«2l/2 e «03i/2\ 



(3) 



It is known since rather long time from the analysis o: 



9 23 and 9 12 are "large" while 13 is "small" (see, e.g., 9 23 = tt/4, 8 12 = sin _1 (va!3) 



the neutrino oscillation data that 



and sin 2 0i3 ^ 0.05 (at 3a). Recently the T2K collaboration reported Q evidence at 2.5a 
for a non-zero value of the "reactor" angle 613. Subsequently the MINOS collaboration also 
reported evidence for a nonzero value of 13 , although with a smaller statistical significance 
4j. A global analysis of the neutrino oscillation data, including the data from the T2K and 
MINOS experiments, performed in [5| showed that actually sin0i 3 7^ at > 3o\ The results 



of this analysis, in which the neutrino mass squared differences Am 2l = Am 2 







2 



and Anijj = Am}, responsible for the solar, and the dominant atmospheric, neutrino oscil- 
lations, were determined as well, are shown in Table [I] The authors of [5] find, in particular, 
the following best fit value and 3a allowed range of sin 2 6* 13 : 

sin 2 9 13 = 0.021 , 0.001 < sin 2 9 13 < 0.044 , (4) 

using the "old" fluxes of reactor v e in the analysis (see jij] for further details). Moreover, it 
was found in the same global analysis that cos 5 = —1 (and sin #13 cos S = —0.14) is clearly 
favored by the data over cos 5 = +1 (and sin # 13 cos 5 = +0.14) [6]. It is interesting to 
note that the best fit value of the Dirac CP violating phase 5, obtained in the analysis of 
the atmospheric neutrino data by the Super Kamiokande (SuperK) Collaboration reported 
in 8], reads 5 = 220° = 1.22-71". A not very different best fit value of 5 was found in the very 
recent global neutrino oscillation data analysis performed in [9]: 5 = — 0.617T (— 0.4l7r), or 
equivalently 5 = 1.397T (1.597r), the value quoted corresponding to neutrino mass spectrum 
with normal (inverted) ordering. It should be added, however, that, except possibly for the 
negative sign of cos 5, the preference of the data for the indicated specific values of the phase 
S are at present statistically insignificant. 

The T2K and MINOS results will be tested in the upcoming reactor neutrino experiments 
Double Chooz jl0|, Daya Bay [11] and RENO [12]. If confirmed, they will have far reaching 



implications for the program of future research in neutrino physics (see, e.g., [13|). 

The experimentally determined values of the solar and atmospheric neutrino mixing an- 
gles 9 12 and # 23 coincide, or are close, to those predicted in the case of tri-bimaximal mixing 

Table I: The best-fit values and 3a allowed ranges of the 3-neutrino oscillation parameters, obtained 
from a global fit of the current neutrino oscillation data, including the recent T2K and MINOS 
results (from [5fl). 
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n 



sin 2 fl 2 T 3 BM 



1/2, sin^™ M = l/3, sin 6 



TBM 

13 



(5) 



with the TBM mixing matrix given by: 



/ 



U 



TBM 



(6) 



\-^/T/6 ) 

It is appealing to assume, following jjjj], that in leading approximation we have (up to 
diagonal phase matrices) C/pmns — ^tbmj the reason being that the specific form of £/tbm 
can be understood on the basis of symmetry considerations. However, in order to account 
for the current neutrino mixing data, and more specifically, for the fact that #13 7^ 0, Utbm 
has to be "corrected". Such a correction can naturally arise given the fact that, as is well 
known, the PMNS matrix receives, in general, contributions from the diagonalisation of the 
neutrino and charged lepton mass matrices, U u and llj, respectively: 



U = U\U V =W$U„Q, 



(?) 



where [15| U\ and U v are V-like 3x3 unitary matrices (see eq. (j2J)) each containing, in 
general, one Dirac-like CP violation phase, $ is a diagonal phase matrix containing, in 



general, two CP violation phases, and Q is define in eq. 
model of flavour we are going to consider in what follows [16 



. In the SUSY 517(5) x T 
17| . U u coincides with Utbm, 
by U}. This 
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- |20| as well 



U v = Utbm, and the requisite "correction" leading to #13 ^ is provided 
possibility was investigated phenomenologically by many authors, see, e.g. 

as Q, ail- 
it has been realized a rather long time ago that the TBM matrix can arise from an 



underlying A4 flavour symmetry 



22j. Nevertheless, 



generating the quark mixing utilizing the 



23|. As a consequence, the incorporation of 



A 4 symmetry encounters a number of difficulties 
the A4 symmetry in GUTs is not straightforward, being plagued with complications. On the 
other hand, the group T' 
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24j . which is the double covering of A4, can successfully account 



for the quark masses and mixing. In |l6l . Il7] a viable SUSY model of flavour, including the 
CP violation in the quark sector, based on the SU(5) x T' symmetry was constructed. The 
model includes three right-handed neutrino fields which possess a Majorana mass term. 
The light neutrino masses are generated by the type I see-saw mechanism and are naturally 



small. The light and heavy neutrinos are Majorana particles. The model is free of discrete 



gauge anomalies 
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26]. In addition to giving rise to realistic masses and mixing patterns 



for the leptons and quarks, the SU(5) x X" model under discussion exhibits a number of 
unique features. In particular, the CP violation, predicted by the model, can entirely be 
geometrical in origin 17). This interesting aspect of the SU(5) x T' model we will consider 



is a consequence of one of the special properties of the group T', namely, that its group 
theoretical Clebsch-Gordon (CG) coefficients are intrinsically complex [27j . More specifically, 
the only source of CP violation, e.g., in the lepton sector of the model is the Dirac phase S. 
The Majorana phases CH21 and 0:31 are predicted (to leading order) to have CP conserving 
values. The Dirac phase is induced effectively by the complex CG coefficients of the group 
T'. The model allows for a successful leptogenesis, the Dirac phase S providing the requisite 
CP violation for the generation of the observed baryon asymmetry of the Universe [28] . As 
a consequence, there is a strong connection between the CP violation in neutrino oscillations 
and the matter-antimatter asymmetry of the Universe. 

In the present article we investigate some aspects of the low energy phenomenology of the 



SU(5) x X" model of flavour 
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17| . which allows to describe in a unified way the masses 



and the mixing of both the quarks and the leptons, neutrinos included. We concentrate first 
on the predictions of the model for the absolute scale of neutrino masses, the neutrino mass 
spectrum and the effective Majorana mass in neutrinoless double beta ((/3/3) ^-) decay. We 
derive next detailed predictions for the rates of the lepton flavour violating (LFV) charged 
lepton radiative decays \i— > e + j, t — >■ e + 7 and r — > jj, + 7. 

The paper is organized as follows. Section 2 is devoted to a rather detailed description of 
the SU(5) x T' model of interest. In Section 3 the predictions of the Model for the neutrino 
masses, the type of neutrino mass spectrum and the (/3/3)oi,-decay effective Majorana mass 
are obtained. In Section 4 we present results on the rates of the LFV charged lepton radiative 
decays \x — >■ e + 7, r— >• e + 7 and r — > // + 7, calculated assuming the mSUGRA scenario of 
soft SUSY breaking. Section 5 contains a Summary of the results obtained in the present 
work. 
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Table II: Field content of the SU(5) x T' model. The three generations of matter fields in 10 and 
5 of SU(5) are in the T3, T a {a = 1, 2) and F multiplets. The Higgses that are needed to generate 
577(5) invariant Yukawa interactions are H§, H— and A45. The flavon fields (j) through Q' are those 
that give rise to the charged fermion mass matrices, while £, 77, and S are the ones that generate 
neutrino masses. The Z\% and Z' 12 charges are given in terms of the parameter u = e"7 6 . 

II. THE SU(5) x T' MODEL 

The content of the chiral superfields of the model considered (including the three gener- 
ations of matter fields, the SU(5) Higgses in the Yukawa sector, and flavon fields) as well as 
their quantum numbers with respect to 577(5), T', and Z12 x Z[ 2 symmetry groups are given 
in Table [TXJ The model includes three right-handed (RH) neutrino fields Nm, I = e, fi, r, 
which are S77(5) singlets, but are assumed to form a triplet of T' . This particle content 
leads to the following Yukawa superpotential up to mass dimension seven: 



(8) 



where 



W TT = y t H 5 T 3 T 3 + —H 5 



1 



+ —y u H 5 T a T b <jj' 3 



W„ = XiNNS + — 
A 6 



y ts T 3 T a ilj( + y c T a T b 

J.Y" 

y s A 45 FT a (PijC + y d H^FT a( p 2 ij' 
H 5 FN(('(\ 2 Z + \ 3 r) 



The UV completion of these operators is discussed in Ref. 



(9) 
(10) 

(11) 



29j. Here the parameter A is the 



scale above which the T" symmetry is exact. The vacuum expectation values of the flavon 
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fields are given by: 



(?) 



M 
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V 1 / 



&A , (4,') 
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, {</>) 







6 A 



(12) 



^ A , m 



(13) 
(14) 



1 

o / vi 

(C) = CoA , (CO = CoA , (v) = VoA , (S) = s A . 

Note that all vacuum expectation values are assumed to be real and they don't contribute 
to CP violation. The Yukawa couplings are also real. The reality of the Yukawa coupling 
constants is ensured by the presence of sufficient number of the complex flavon fields which 
allows to absorb the complex phases in the Yukawa coupling constants by phase redefinitions 
of the fields. 

The superpotential gives rise to the following mass matrix for the up-type quarks, 

\ 



o 







-^jyo 0o + (1 - |)0o 



V 



VtV u , 



2/VoCo 1 

and the following down-type quark and charged lepton mass matrices, 



(15) 



/ 



M d 



\ 



\ 



ydVd<Po 



(l + ^)0o^o 

-(i-i)M' Q VoCo o 

00^0 0o^o Co J 

-(1 - i)(j> Q % 00^0 

1 + i)M' -3^oC 0o^o 
Co J 

The model describes successfully the quark masses and mixing, as discussed in detail in 



(16) 



ydVd<Po 



(17) 



16J. 



The presence of the complex elements in M u and Md allows to account for the observed 



CP violation in the quark sector as well 
sector. 



30|. We concentrate in what follows on the lepton 
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In the lepton sector, the superpotential leads to the following Dirac neutrino mass matrix, 



M, 



D 



2Co + ?7o -Co -Co 
-Co 2^o -Co + Vo 
\ -Co -Co + ^70 2£ J 

and to the following RH neutrino Majorana mass matrix 



CoCqVu = Y v v u 



(18) 



^10 0^ 



M 



RR 



1 

1 



s A. 



(19) 



We note that the complex CG coefficients appear in the product rules involving the 
spinorial representations 2, 2', 2" of T'. Because (T 1 ,T 2 ) transform as the spinorial rep- 
resentation 2, the charged fermion mass matrices, M u , Md, M e , are complex. On the 
other hand, the neutrinos involve only the vectorial-like representations 1, 1', 1", 3 of T'. 
Therefore the neutrino Dirac and Majorana mass matrices are real and thus CP conserving. 

Note also that the Dirac neutrino mass matrix, Md, is real and symmetric. As can be 
easily shown, it is diagonalized by the TBM matrix, 



^tbm m d U TBM = M D iag = diag(3£, + rjo, r} , 3£ - Vo) 



(20) 



where all elements in the diagonal matrix M D iag are real. 

The RH neutrino Majorana mass matrix Mrr is diagonalised by the unitary matrix S: 



S T MrrS = D N = diag(M 1 ,M 2 ,M 3 ) = s Adiag(l, 1,1), Mj > 0, j = 1,2,3, (21) 



where 



/ 



s 



\ 



1 

1/V2 -i/V2 

^0 l/y/2 l/V2 J 

and Mj are the masses of the heavy Majorana neutrinos Nj (possessing definite masses 



(22) 



S^Nm + Sj t C(N jR y = C(NjY , j = 1, 2, 3 



(23) 



C being the charge conjugation matrix. Thus, to leading order, the masses of the three 
heavy Majorana neutrinos Nj coincide, Mj = SqA = M, j = 1, 2, 3. It follows from eq. (12 ip 
that S'*S'^ is a real matrix, so S*S^ = SS T . 
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The effective Majorana mass matrix of the left-handed (LH) flavour neutrinos, M v , which 
is generated by the see-saw mechanism, 

M v = -M D M^ l R M T D , (24) 

is also diagonalized by the TBM matrix, 

Utbm M uUtbm = diag((3£ + Vo) 2 , Vo, -{~Ho + Vo) 2 ) ° g ° A " = Q diag(m h m 2 , m 3 ) Q T . 

(25) 

Here Q is a diagonal phase matrix, Q = diag(l,l,±i), and rrik > 0, k = 1,2,3, are the 

masses of the three light Majorana neutrinos: 

(X + 3Z) 2 X 2 (X-3Z) 2 , x 

mi " M ' m2 "M' m3 " M ' (26) 

where X = 77o(CoCo' y «) an d ^ = ^(CoCo^u)- ^ follows from eqs. (j7|) and ( 123]) that U U Q = 
Utbm Q- Thus, we have Q = Q and the Majorana phases a 2 i and 031 are predicted (to 
leading order) to have CP conserving values: a 2 \ = and 031 = ±7T. 

One special property of M v is that it is form diagonalizable |30|. In other words, regardless 
of the values of £ and i] , M v is always diagonalized by the TBM matrix, U TBM . 

The charged lepton mass matrix M e , eq. (TTTj) . is diagonalised, in general, by the by- 
unitary transformation: M e = U e M^V^ R , where V e R and U e are unitary matrices and M e d = 
diagim ei m iJil m T \ mi being the mass of the charged lepton I, I = e,fi,r. Thus, the matrix 
U e , which enters into the expression for the PMNS matrix, U = U^U P , digonalises the matrix 
M e Mh Ul(M e M\)U e = (M e d ) 2 . The 3 x 3 unitary matrix U e of interest can be parametrised, 
in general, as U e = $ V e (9\ 2 , 6f 3 , 0f 3 , S e )P, where V e (8^ 2 , 8f 3 , 8 23 , S e ) has the same form as the 
matrix V in eq. (j2J) (with 6^- and S replaced by 9L and 5 e ) and $ and P are diagonal phase 
matrices containing two and three phases, respectively. The phases in P are unphysical and 
we are not going to consider them further. 

Due to the SU(5) symmetry, the matrices and M e depend on the same three real 
parameters and two phases, the latter having the values ±7r/4. The three real parameters 
can be fixed by requiring that and M e reproduce correctly the down-quark and charged 
lepton mass ratios rrid/m s , m e jm^ etc., the Cabibbo angle 9 C and other quark mixing 
observables. The model accounts successfully for the quark masses and mixing and for the 
charged lepton masses flf] . \\\ . In particular, the well-known relations sin# c = ^J~m^rn Sl 
m e / m [i — ra<f/(9m s ), etc. are fulfilled. 



Fitting the indicated quark sector observables and charged lepton masses one finds that 



17j | two of the three angles in the matrix U e are extremely small, 



sin ^ 1.3 x 10" 



sin 9 e 23 ^ 1.5 x 10" 



(27) 



while the third satisfies: 



sm6 e 



12 



sin# r 



(28) 



It follows from the quoted results that, to a very good approximation, we can set 9\ 3 = 
#23 = in the expression for U e . In this approximation we have U e = $i?i2(0f 2 )> where 
$ = diag(l, e iif , 1) and 



/ 



R 



12^12 J 



cos 9% sin^fo 



\ 



sin 9\ 2 cos 9f 2 



V 



(29) 



1 J 

Comparing the expressions in the left-hand and right-hand sides of the equation M e MJ = 
U e {M^) 2 Ul and assuming, without loss of generality, that ^oV-'o^oCo > and cos9f 2 sin#f 2 > 
0, we find that 



7T 

tp= T 

In the approximation we are using the PMNS matrix is given by: 



jjSU(5)xT' = jjl 



R 



12(^12 



§*{v)UtbmQ 



( v/273cf 2 + y/TJEs^e-^ y/T/3(c e 12 - s\ 2 e~^) y/lfisf^ \ 



^2/3st 2 -y/T/6c e 12 e-^ y/T/3(s e 12 + c e 12 e-^) -^lj2c\ 2 e 



■up 



\ 



■y/T/6 



(30) 



(31) 



Q ,(32) 



where c f 9 = cos#i 2 , s i2 = sin0* 2 . It follows from the above expression and eqs. ([I]) - ([3]) that 
15j , Il9l . |20| up to corrections of the order of sin 2 0^ 2 , #23 takes its TBM value, #23 = 7r/4, 
and that to leading order in sin 9\ 2 we have 



- * sin 



and 



12 1 2y/2 
sm 6»i2 = - - - sin 9 12 cos ip = — sin 9 13 cos ip . 



Using eqs. ([28]) and ([30]) we get (l7|: 



sin 6»i3 = — y= sin 9 C , 



(33) 



(34) 



(35) 
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and 



.2, ~ 1 V2 



sm 2 6 12 = --^r sin^ c . (36) 



Thus, in the model considered the CHOOZ angle #13 is predicted to be rather small: using 
sin# c = 0.22 we get from eq. ( 1351) . sin 613 = 0.052. From a numerical analysis in which the 
higher order corrections were also included one finds jl7| sin #13 = 0.058. This value lies 
in the 3a interval of allowed values of sin# 13 , determined in the global analysis of the 
neutrino oscillation data. The correction to the TBM value of sin 2 9i 2 given in eq. (136|) . is 
negative. The value of sin 2 9± 2 = 0.299 predicted by the model, lies within the la allowed 
range, found in the global data analysis 

It is possible to relate also the Dirac CP violating phase 5, present in U, eqs. (CEJ) - ([3]) , with 
the phase (p in U' . It proves convenient first to multiply the elements of the first row and of 
the first and of the second columns of the PMNS matrix in eq. (132|) by (-1), and the elements 
of the second row by (— e %(p ). Multiplying by (-1) {{-e %Lp )) the first (second) row is equivalent 
of redefining the phase of the electron (muon) field in the weak charged current. Changing 
the signs of the elements of the first and of the second columns of the PMNS matrix in eq. 
( 1321) can be compensated by multiplying the matrix Q containing the two Majorana phases 
by diag(— 1, — 1, 1): diag(— 1, — 1, 1)Q = diag(— 1, — 1, ±z) = (—l)diag(l,l,^fi) = (—1)Q*, 
where the overall factor (-1) in the matrix Q* has no physical significance and will be dropped 
in our further discussions. After this simple manipulations the matrix in eq. (132]) takes a 
form which is similar to that of the standard parametrisation of the PMNS matrix, eqs. (QQ) 



U' 



Q* , (37) 



where 



v/ = j±tt. (38) 

We note that the phases of the U' el and U' e2 elements of U', K e % and n e2 , are exceedingly 
small: K el = s e 12 sinp'/2 = (-0.0259) = (-1.5°), K e2 ^ -s^siny?' ^ 0.0519 = (3.0°). Thus, 
the imaginary parts of U' el and U' e2 are much smaller than their real parts and we have 
U' ek = \U> k lk = 1,2. 
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Comparing the real and imaginary parts of the quantity U^UjjiUesU^, calculated using 



eqs. (pp) - (j3J), with those obtained utilizing eq. (137]) (see, e.g., |20j) and assuming that the 
Dirac phase 6 lies in the "standard" interval [0, 2ir} we find 

5 

5 = if' = if + 7T ^ -n. (39) 

Note that the sign of cos 5 and the value of 5 are compatible with those suggested by the 
current neutrino oscillation data. Substituting tp with (5— 71") in eq. ( l36l) and using the results 
obtained in [20] on the values of cos 5, allowed by the existing data on sin 2 9u and sin #13, we 
get for sin #13 = 0.058, predicted by the model, and the 3a (2a) interval of experimentally 
allowed values of sin 2 #12 : 

- 1 < cos 5 £ 0.4 (0.1) . (40) 

Thus, in the model considered, 5 = and, more generally, the values of cos 5 from the 
interval 0.4 < cos 5 < 1, are excluded at 3a. 

As is well known, the quantity Jcp = Im^^U^iU^U*^ is the rephasing invariant associ- 
ated with the Dirac CP violating phase 5 in the PMNS matrix. It determines the magnitude 
of CP violation effects in neutrino oscillations [32] and is analogous to the rephasing invariant 
associated with the Dirac phase in the Cabibbo-Kobayashi-Maskawa quark mixing matrix, 



introduced in 



331 ] . In the model considered the rephasing invariant Jqp is given to leading 



order in sin 9\ 2 by 



J cp = - sin^f, sin 5 = — sin# c sin 5 = = sin# c . (41) 

° 6 12 18 18^2 



Finally, we give the expression for the PMNS matrix U' obtained numerically in 



13], in 



which the higher order correction in sin#^ 2 , as well as the corrections due to the nonzero 



values of sin^g ~ 10 5 and sin 6^3 ~ 10 4 , are all taken into account: 



jjSU(5)xT' _ jji 



\ 



Q. (42) 



0.838e iL4 ° 0.543e i6 - 5 ° 0.058e" i22L80 
- 0.363e- i3 - 53 ° 0.610e i7 2 ° 0.705e i4 ° 
y 0.408 -0.577 0.707 j 

It follows from this expression of U' that the precise values of sm8 13 , 5 and Jqp, predicted 
by the model, read: sin# 13 = 0.058, 5 = 221. 80 and Jcp = —9.66 x 10~ 3 . The predictions 
of the SU(5) x T' model of interest for sin ^13, sin 2 612, 8 and Jqp can be tested directly in 
the upcoming neutrino oscillation experiments. 
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III. NEUTRINO MASSES, SPECTRUM AND THE (/3/3) 0i ,-DECAY EFFECTIVE 
MAJORANA MASS 



In the model considered we have cos 26 12 > and therefore Am^ should also be positive. 
It is easy to show then that we also have Am\ = Am^ > 0, i.e., that in this model 
only light neutrino mass spectrum with normal ordering (or normal hierarchy) is possible. 
Indeed, using the expressions given in (1261) one finds that there do not exist values of X 
and Z for which the condition which defines the spectrum with inverted ordering, namely, 
Am\ = Amgj < 0, or m 2 > mi > m 3 , is satisfied. 

The neutrino masses in the model considered, as it follows from eq. (I2"6"j) . depend on two 
parameters: X/yM and Z/yfM. The latter can be determined, e.g., using the data on 
Am\ and Am^. Thus, the absolute values of the three light neutrino masses m 12)3 in the 
SU(5) x T' model under investigation are fixed by the values of Am^ and Am^. 

It proves convenient to use Am^ and the ratio r = Am^/Am^, instead of Am^, to 
determine the masses mi )2i3 . We have for the 3a allowed range of the indicated ratio: 

Ami Am 2 ,, X 4 - (X + 3Z) 4 , . 

r = ir^r = iHr = tt? ^ — tt? tzt* = 0-032 ± 0.006, 43 

Am\ Am\ x (X - 3Z) A - (X + 3Z) 4 ' v ; 

Using the fact that, e.g., ml = Am\ x + m\ and eq. ( 1261) . we can express the heavy Majorana 
neutrino mass M in terms of Am^, X and Z: M 2 = (X A — (X + 3Z) 4: ) / (Amli) . Substituting 
this expression for M in eq. (12611 . we get for mf: 



2 ^ m ll 2 Am 21 2 Am 21 

X \ 4 /X + 3Z\ 4 / X \ 4 /X + 3Z x4 



X + 3ZJ \ X J \X-3ZJ \X-3Z 

'(44) 

Given r, eq. (l4"3"|) implies a relation between the parameters X and Z. As can be shown, to 
first order in r, this relation reads: 3Z ~ (20r — 2)X. A numerical analysis performed by 
us showed that the values of the light neutrino masses, calculated using eqs. ( |26|) and (JJ3J) 



and the values of Am^ and r as input, are reproduced with a remarkable precision when 
calculated employing the approximate relation between Z and X given above instead of the 
exact one implied by eq. (|43p . Utilizing the relation 3Z ~ (20r — 2)X and eq. (|44p allows 
to express m\ 2 3 as simple functions of Am| x and r: 

Ami, (20r - l) 4 , Am2, 2 A£ (20r - 3) 4 , . 

2 ^ 21V ; 2 ^ 21 2 ^ 21V ; 45 

1 i_(20r-l) 4 ' 2 l-(20r-l) 4 ' 3 1 - (20r - l) 4 1 ; 
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We also have: 

M \jbm\ x = X 2 {\- (20r - l) 4 ) * = X 2 , (46) 

where we have neglected (20r — l) 4 = 0.017. For, e.g., M = 10 12 GeV, using X = ^(CoCo^m); 
v u = 174 GeV and Am 2 21 = 7.59 x 10~ 5 eV 2 we obtain: (r/oCoCo) 2 - 2 -89 x 10~ 4 . Taking into 
account that Z = £o(CoOX)> we find also £o = (20r - 2)r/ /3 = -0.453r/ . For M = 10 12 
GeV we have: &C0C0 - -7-74 x 10~ 3 . 

In Fig. [1] the light neutrino masses 777.1,2,3 > are plotted as functions of the ratio r, 
with Am^i set to its best fit value. As Fig. [1] indicates, in the model with approximate 
SU(5) x T' symmetry under study, the light neutrino masses are allowed to vary (due to 
the uncertainties in the experimentally determined values of Am^ and r) in rather narrow 
intervals around the values mi = 1.14 x 10~ 3 eV, m 2 = 8.78 x 10~ 3 el^, m 3 = 4.89 x 10~ 2 eV. 
This, implies, in particular that the model provides specific predictions for the sum of the 
light neutrino masses as well as well for the effective Majorana mass in neutrinoless double 
beta ((/3/3)o I y-decay), |<m>| . 

For the sum of the three light neutrino masses we obtain: 

(Am 2 \ 5 
l-(20r-l)* j [1 + (20r - l) 2 + (20r - 3) 2 ] (47) 

Numerically we get using the best fit values of Am^ and r: 

mi + m 2 + m 3 ^ 5.9 x 10~ 2 eV , (48) 



To leading order in sin 9f 2 , the (/3/3) j,-decay effective Majorana mass (see, e.g., 34]), 
|<m>| = I J2l=i( U PMNs)tk m k\, is § iven b y : 



~ (2mi + m 2 ) - ~ (m 2 - m x ) sin fl^e"* + i m 3 (sin #y 2 e^"*) 



(49) 

where we have employed the expression for Upmns given in eq. (132|) . Using the numerical 
values of m 12i 3 corresponding to r = 0.032 quoted above, eqs. (128]) and ( |30|) and the value 
of sin# c = 0.22, we find that the term ex isimp and that cx m 3 in |<m>| give negligible 
contributions and that 

|<m>| S 3.4 x 10~ 3 eV. (50) 
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Figure 1: The light neutrino masses mi (solid line), mi (dashed line) and m.3 (long-dashed line) 
as functions of the parameter r for Am^ = 7.58 x 1CP 5 eV 2 . The three vertical lines correspond, 
from left to right, to r = 0.026, 0.032, 0.038, i.e., to the 3a minimal, best fit and 3a maximal values 
of r. For r = 0.032 we have: mi = 1.14 x 10" 3 eV, m 2 = 8.788 x 10" 3 eV, m 3 = 4.89 x 10~ 2 eV. 



IV. THE LFV DECAYS L 



+ 7 



A. The Decays 



ij + 7 in SUSY theories 



In the minimal extended Standard Theory with heavy Majorana right-handed neutrinos, 
the simultaneous presence of neutrino and lepton Yukawa couplings, Y v and Y e , leads to 
lepton flavour violation. In this scenario, the LFV decay rates and cross sections are strongly 
suppressed by the ratio 
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360 l(E,=2 3 U llj U eA m h)/M^\ 2 , M w 80 GeV being the W= 



boson mass, leading to, e.g., BR(/i — » e7)< 10 -47 . This renders the LFV decays and 
reactions unobservable in the ongoing and future planned experiments. The situation is the 
same in the non-SUSY type I see-saw models in which the heavy Majorana neutrino masses 
Mfc are only by few to several orders of magnitude smaller than the GUT scale Mx ~ 2 x 10 16 



GeV 



37 



38]. 



The present experimental upper bounds on the rates of the LFV decays 



m £i > m e i x = e, £ 2 = n, £3 = t, are given by 



39 



40 1 



+ 7, 



BR(fi -> erf) < 2.4 x 10" 12 , BR(p -)■ 3e) < 1 x 10" 12 , BR(t ->■ // 7 ) < 4.4 x 10" 8 . (51) 

The first was obtained recently in the MEG experiment at PSI and is an improvement by 
a factor of 5 of the previous best upper limit of the MEGA experiment, published in 1999 
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Figure 2: Feynman diagrams of the dominant contribution to the — > £j + 7 radiative decay 
amplitude in the mass insertion approximation. \A corresponds to charginos or neutralinos and Lj 
are the slepton doublets. The photon can be emitted from the chargino or the slepton lines. 



41]. The projected sensitivity of the MEG experiment is 39]: 



BR(p ->■ erf) ~ KT 13 . (52) 



In the supersymmetric (SUSY) theories these decay rates can be largely enhanced due to 
contributions from the slepton part of the soft SUSY breaking Lagrangian, C so f t : 

- C sof t = {m\) i3 l\Lj + {ml) i3 e Ri e Rj + (ml) i3 P Ri u Rj 

+ {{A e ) t3 H d e R ~L 3 + (A u ) t3 H u v* m L 3 + h.c.^j , (53) 

where m| and m| are the left-handed (LH) and right-handed (RH) charged slepton mass 
matrices, respectively, m? is the right-handed sneutrino soft mass term, A e and A v are 
trilinear couplings and and H u are the two Higgs doublet fields present in the SUSY 
theories. Non-zero off-diagonal elements in the slepton mass matrix could induce lepton 
flavour violation, but only relatively small values could satisfy the upper bounds quoted 
above. Soft-breaking terms are, however, subjected to renormalization through Yukawa and 
gauge interactions in such a way that LFV is induced in the slepton mass matrix at "low 
energies". Indeed, in addition to a LF conserving part, the renormalisation group (RG) 
equation for the left handed slepton mass matrix present, in general, off-diagonal terms 
which are a source of LFV. 

The indicated generic possibility is realized in the SUSY (GUT) theories with see-saw 



mechanism of neutrino mass generation 



42|. If the SUSY breaking occurs via soft terms 



with universal boundary conditions at a scale Mx above the RH Majorana neutrino mass 
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scale Mr, Mx > Mr, as in the so-called minimal supergravity (mSUGRA) scenario 43], 
the renormalisation group effects transmit the LFV from the neutrino mixing at Mx to the 
effective mass terms of the scalar leptons at Mr, generating new LFV corrections to the 
flavour-diagonal mass terms. For slepton masses of a few hundred GeV, the LFV mass cor- 
rections at Mr are typically of the order of a few GeV and thus are much larger than the light 
neutrino masses rrij . As a consequence (and in contrast to the non-supersymmetric case), 
the LFV scalar lepton mixing at Mr generates additional contributions to the amplitudes 
of the LFV decays and reactions which are not suppressed by the small values of neutrino 
masses. As a result, the LFV processes can proceed with rates and cross sections which are 
within the sensitivity of presently operating and future planned experiments 



also, e.g., 



42 



44] (see 



45M49| and the references quoted therein). 



In the following discussion we will assume the commonly employed mSUGRA SUSY 



breaking scenario 



43] . In this scenario the flavour is assumed to be exactly conserved at 



the GUT scale, Mx ~ 2 x 10 GeV, by the soft SUSY breaking terms. More specifically, 
it is assumed that at the scale Mx the slepton mass matrices are diagonal in flavour and 
universal, the trilinear couplings are proportional to the neutrino and charged lepton Yukawa 
couplings Y v and Y e , respectively, and the gaugino masses have a common value: 



[m 



Mb 



2\ _ r 2 
ij ~ °ij m 



Y u a m , A e = Y e a m 



My 



w 



Mg = mi/2- 



(54) 



Hence the parameter space of interest of mSUGRA is determined by: 



mi/2, mo, A , tg/3. 



(55) 



tg(3 being the ration of the vacuum expectation values of the Higgs fields H u and 

In the leading-log approximation, the branching ratios of the LFV processes i± — > i j + 7 



(m^ > mi ) is given by 4J| (see also, e.g. 



45 



47 



B(£ t t s + 7) = ^f^p- B (m , m 1/2 ) 



^ M x M x f 

k 



tan 2 /5. (56) 



Here Gf is the Fermi constant, a em ~ 1/137 is the fine structure constant, is the mass of 
the heavy Majorana neutrino A^, Y v is the matrix of neutrino Yukawa couplings in the basis 
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m [GeV] m 1/2 [GeV] 



Figure 3: The scaling factor Bo(mo, m^) as function of uiq (left panel) and m\/2 (right panel) for 
Aq = 7mQ. Left panel: the values of used are 200 GeV (solid line), 400 GeV (long-dashed 

line), 600 GeV (medium long-dashed line), 800 GeV (short-dashed line). Right panel: the values of 
mo used are 50 GeV (solid line), 100 GeV (long-dashed line), 200 GeV (medium long-dashed line), 
400 GeV (short-dashed line). 

in which the Major ana mass matrix of the RH neutrinos and the matrix of charged lepton 
Yukawa couplings are diagonal, and J40| T(/j — > evv)/T tot (n) « 1, T(r — > euP)/T tot (r) ~ 
0.1785, T(r — » [ivi)) /T tot (r) m 0.1736. The scaling function B Q (m Q ,mi/2) contains the 
dependence on the SUSY breaking parameters: 



B Q (m ,mi/ 2 ) 



I (3 + al)m 



G 2 F m% 



(57) 



8tt 2 

The effective SUSY mass parameter mj that appears in equation (156"]) can be approximated 



by 



47J: 



m 



s ~ 0.5mX/ 2 K + 0-6mf /2 ) 2 (58) 



This analytic expression was shown to reproduce the exact RG results for B(&i — >■ ^+7) with 
high precision. In Fig. [3] we illustrate the dependence of the scaling function B (m Q ,mi/2) 
on the parameter itlq {m-y/2) for four values of mi/j (mo). 

B. Predictions of the SU{5) x T Model 

As we have seen, in the SU(5) x T' model of flavour we are considering the three heavy 
Majorana neutrinos Nj are, to leading order, degenerate in mass: Mj = SoA = M > 0, 
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j = 1,2,3. The higher order corrections to the masses Mj lead to exceedingly small effects 
in the £i — > £j + 7 decay rates and we will neglect them. In this case the branching ratios 
B(£i — > £j + 7) of interest, as it follows from eq. ( |56|) . depend on the quantity | (Y^l^t)^- 1 2 . 

As can be shown, in the model under study, the matrix of neutrino Yukawa couplings Y v 
entering into the expressions for B[£i — > £j + 7) is related to the Dirac neutrino mass matrix 
Mo, eq. (TT8]) . as follows: 

Y u = — U\ M D S (59) 

Vu 

where U e is the matrix diagonalising M e Mj", M e being the charged lepton mass matrix, and 
5* is determined in eqs. ([211) and (1221) . In the basis in which the RH neutrino Majorana 
mass matrix and the matrix of charged lepton Yukawa couplings are diagonal, the Majorana 
mass term for the LH flavour neutrinos, generated by the see-saw mechanism, is given by: 

M u = -v 2 u Y U D N X T = U' D u (U'f , (60) 

where D u = diag(mi, 1712, WI3) and 

U' = UlU TBM Q. (61) 

Equation (1601) . as is well known, allows to express Y v in terms of U', D„, Dn and a complex, 
in general, orthogonal matrix 



45] R, R R = RR = 1: 
Y v = — U'y^Ry^. (62) 

V u 

From eqs. (159]) - (162]) and ( ]20]) . we obtain the following expression for the matrix R: 



R = v u (V^)' 1 Q* K ag U% BM S (^D^r 1 . (63) 

Using the explicit form of Q = diag(l, 1, ±i) and of Utbm and S, eqs. ([6]) and (122]) . it is 
not difficult to show that R is a real matrix: R* = R. Taking into account this result and 
recalling that = M diag(l, 1, 1), we get for the quantity of interest Y u Yj: 

Y V Y} = ^U'D V {U')\ (64) 

V u 

where the PMNS matrix U' = U^UtbmQ, see eqs. (J3TJ), ([32]), (JHZJ) and @2}. Thus, the 
£% £j + 7 decay branching ratios in the SU (5) x T' model under investigation depend 
on the mass of the heavy Majorana neutrinos through the factor (M(lnMx/Mfc)/t; u ) 2 , do 
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not depend on the matrix R, and their ratios are entirely determined by the elements of the 

PMNS matrix and the light neutrino masses. More specifically, using the general expressions 

I I 

for the elements of the PMNS matrix and the unitarity of the latter we get [48]: 

K^yJ)^) = [U^rrijU*^ ~ |A 2 i si 2 c 23 ci 2 + A 31 s 23 si 3 e tS \ , (65) 

|(Y^ Y"J) re | = | U Tj rrijU*^ ~ |- A 21 s 12 ci 2 s 2 3 + A 31 C23Si 3 e I<5 | , (66) 

\{Y v Y}) Ttx \ = \U Tjmj U^\ ~ |A 3 i c 23 s 23 | , (67) 

where 

M R Ami , x 

Ay = -4 * J . 68 

V* TUi + Ulj 

In eqs. (165]) - (IBTj) we have used the approximation c 13 = 1 and have neglected c 2 2 A 21 , 
s 2 2 A 21 and higher order terms oc A 21 with respect to A 31 . 

We define the "double" ratios of the branching ratios B(£i — > £j + 7) as follows: 

(69) 

The double ratios of interest are given by: 



\%v2U 2 D/O1/On ,j(^j) 



|2 



M21/31) g lv - V^' , #(21/32) S lv - V . (70) 

|(nrj) re | 2 ' IC^WI 2 

C. Numerical Results 

In Figs. H]- [7] we present results on the branching ratios BR(fi — > e + 7), BR(t — )■ e + 7) 
and BR(r — > /i + 7), predicted by the SU(5) x T" model. The branching ratios are calculated 
using eq. flS>D with M fc and M x set to M fc = M = 10 12 GeV and M x = 2 x 10 16 GeV. The 
values of the mSUGRA parameters tan (3 and A are chosen from the intervals tan = 3-i-50, 



50 



51] and 



Aq = 0-i-7mo, which are compatible with the constaints obtained by the ATLAS 
CMS 52] experiments at LHC. The values of the other two relevant mSUGRA parameters, 



m and mi/ 2 , are chosen from intervals favored by the global data analysis performed 53]. 
The data set used in this analysis includes in addition to the results of the ATLAS and the 
CMS experiments, the data on the muon [g — 2), on the precision electroweak observables, 
on 5-physics observables, astrophysical data on the cold dark matter density, as well as the 
limits from the direct searches for Higgs boson and sparticles at LEP. Based on the results 
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obtained in 53], mo and 7711/2 are chosen from, or to vary in, the intervals 50 GeV < mo < 400 
GeV and 300 GeV < mi/2 < 800 GeV. For the values of the lepton mixing angles and the 
Dirac phase 5 we use those predicted (to leading order) by the 577(5) x T' model (unless 
otherwise stated): sin 2 #23 = 0.5, sin 2 12 = 0.299, sin #13 = 0.058, 5 = 57r/4. The neutrino 
masses used as input are obtained using the best fit values of AmL = 7.58 x 10 -5 eV 2 and 



0.032. 
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Figure 4: The branching ratios BR(^ — > e + 7) (left), BR(r — > e + 7) (center) and BR(r —7-^ + 7) 
(right) as function of m for Am| = 7.58 x 10~ 5 eV 2 , r = 0.032, 6 = 5vr/4, sin(9i3= 0.058, 
tan/3 = 28, Aq— 7mo, and three values of m 1( / 2 : 400 GeV (solid line), 600 GeV (dashed line), 800 
GeV (dotted line). The horizontal line corresponds to the MEG bound BR(/i — > e + 7)< 2.4 x 10~ 12 . 
See text for details. 

As the Figs. H]- [7] show, for the values of M = 10 12 GeV and tan (3 = 28 considered and 
mo from the interval (50 - 400) GeV, the predictions for the branching ratios BR(fi — > ey), 
BR(r — > erf) and BR(t — > fij) are very sensitive to the values of m^: when the latter 
increases from 300 GeV to 800 GeV, the branching ratios decrease by approximately 2 orders 
of magnitude. For the fixed value of A = 430 GeV, we have BR(fi — > e^) ^ 10~ 13 provided 

mo ^ 300 GeV and ^ 400 GeV (Figs. [6]and[7]). The ranges of values of mo and mi/2, for 
which 10~ 13 ^ BR(fi — > e^) < 2.4 x 10~ 12 , is very different in the case when A oc m . For 
A = 7m ( Figs. IH-E]), for instance, BR(fi — > e*y) satisfies the MEG upper limit and is in 
the range of sensitivity of the MEG experiment for m lying in the interval m = (100 — 300) 
GeV if 600 GeV ^ mi/2 ^ 800 GeV. For the values of the parameters used to obtain Figs. 
H-0 we find that BR(r /1 + 7) £ 10~ 9 and BR(t e + 7) £ 3 x 10 -11 . 
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Figure 5: The branching ratios BR(fj, — > e + 7) (left), BR(t — > e + 7) (center) and BR(t —7-^ + 7) 
(right) as function of m^j for three values of tuq: 50 GeV (solid line), 150 GeV (dashed line), 250 
GeV (dotted line). The values of the other parameters used are the same as those quoted in the 
caption of Fig. [U The horizontal line corresponds to the MEG bound BR(^ — > e + 7)< 2.4 x 10~ 12 . 
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Figure 6: The branching ratios BR(// — > e + 7) (left), BR(r — > e + 7) (center) and BR(r — > /x + 7) 
(right) as function of mo for Aq = 430 GeV and three values of m^: 400 GeV (solid line), 600 
GeV (dashed line), 800 GeV (dotted line). The values of the other parameters used are the same 
as those quoted in the caption of Fig. UJ The horizontal line corresponds to the MEG bound 
BR(n ^ e + 7) < 2.4 x 10~ 12 . 

The jj, — > e + 7 decay branching ratio BR(fi — e + 7) can exhibit very strong dependence 
on the value of the angle # 13 if the Dirac phase 5 = tt, and on the Dirac phase 5 in the case 
of sin #13 = 0.058. This is illustrated in Fig. [8] and can be easily understood on the basis 
of the analytic expression for | (V^Vj ) /ue | given in eq. (165|) . As Fig. [8] shows, for the values 
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Figure 7: The branching ratios BR(/i — > e + 7) (left), BR(r — > e + 7) (center) and BR(r —7-^ + 7) 
(right) as function of mi/2 for Aq = 430 GeV and mo = 150 GeV (solid line), 250 GeV (dashed 
line). The values of the other parameters used are the same as those quoted in the caption of Fig. 
UJ The horizontal dashed line corresponds to the MEG bound BR(fi — > e + 7) < 2.4 x 10~ 12 . 
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Figure 8: The branching ratio BR(/x — > e + 7) as function of #13 (left panel) and of the Dirac phase 
5 (right panel) for r = 0.032, m = 230 GeV, m 1/2 = 400 GeV, ^ = 430 GeV and tan j3 = 28. In 
the left panel, BR([i — > e + 7) is plotted for three values of 5: tt/2 (solid line) , ir (dashed line), 57r/4 
(dotted line). The horizontal line corresponds to the MEG bound BR(fi — > e + 7) < 2.4 x 10~ 12 . 
The results for BR(/j> — > e + 7) shown in the right panel are obtained for sin #13 = 0.058. The solid 
vertical line in the right panel corresponding to 5 = bit /A. 



of 5 = 5tt/A and sin #13 = 0.058, predicted by the model, BR(fi — > cy) is somewhat smaller 
than the maximal value it can have as a function of S and sin# 13 . 
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We have studied also the dependence of the [L — > e + 7 decay branching ratio BR([i — > 
e + 7) on the parameter r. The results of this study are illustrated in Fig. [9] for three 
values of sin 013, sin 2 # 13 = 3.4 x 10" 3 , 0.01, 0.02, and three values of 5, 5 = vr/2, vr, 5vr/4. 
For illustrative purposes the ratio r is varied in the interval < r < 0.05, which is wider 
than the current 3a range of allowed values of r, 0.026 < r < 0.038. Figure exhibits in a 
different way the sensitivity of BR{^i — > e + 7) to the values of #13 and S. 
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Figure 9: The dependence of the branching ratio BR(/x — > e + 7) on the parameter r for mo = 230 
GeV, mi/2 = 400 GeV, = 430 GeV and = 28 . The «, mlddl e and right pan e ls M e 

obtained for sin 2 &13 = 3.4 x 10 3 [17l |. 0.01 and 0.02, respectively. The values of S used are ir/2 
(solid line), tt (dashed line) and 57r/4 (dotted line). The vertical lines correspond respectively to 
r =0.026, 0.032, 0.038. The horizontal line shows the MEG bound BR(/j, -> e + 7) < 2.4 x 10" 12 . 
See text for details. 

In Fig. [10] we present results the double ratios i?(21/31) and i?(21/32) predicted by the 
model, as a function of the lightest neutrino mass mi. We recall that the double ratios 
i?(21/31) and i?(21/32) in the model considered do not depend on the heavy Majorana 
neutrino mass M and on the mSUGRA parameters: they are entirely determined by the 
values of the neutrino masses, the neutrino mixing angles and the Dirac CP violating phase 
5. All these neutrino parameters have essentially definite values in the SU(5) x X" model 
considered. Thus, the values of the double ratios i?(21/31) and 72(21/32) are predicted by 
the model with relatively small uncertainties, as is also seen in Fig. [lOj we have numerically 

#(21/31)^0.21, #(21/32) =■ 7.4 x 10~ 3 . (71) 

These values are one of the characteristic predictions of the 577(5) x T" model of flavour 
under investigation. It is interesting to note also that in the model with 577(5) xT" symmetry 
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Figure 10: The double ratios R(21/31), left panel, and R(21/32), right panel, as functions of 
the lightest neutrino mass, mi. The three vertical lines in each plot are positioned at mi = 
2.06 x 10~ 3 , 1.14 x 10~ 3 , 5.0 x 10~ 4 eV, and correspond to r = 0.026, 0.032, 0.038 and the best fist 
value of Am^i = 7.58 x 10~ 5 eV 2 . The solid lines are obtained with the best-fit values of Am^ 
and Amgj given in Table HI while the dashed lines indicate the 3a allowed ranges. 

considered, the r — > /i + 7 decay branching ratio BR(r — > /x + 7)/ BR(t — > ev T v e ) can be 
bigger than the \i — > e + 7 decay branching ratio BR(fi — > e + 7) by a factor ~ 10 2 . 



V. CONCLUSIONS 



In the present article we have investigate certain aspects of the low energy lepton phe- 

17| . and 



nomenology of the SUSY SU(5) x T' model of flavour which was developed in [16 
which allows to describe in a unified way the masses and the mixing of the quarks and the 
leptons, neutrinos included, as well as the CP violation both in the quark and lepton sectors. 
The model includes three right-handed neutrino fields which possess a Majorana mass term. 
The light neutrino masses are generated by the type I see-saw mechanism and are naturally 
small. The light and heavy neutrinos are Majorana particles. The heavy Majorana neutri- 
nos Nk are predicted (to leading order) to be degenerate in mass: = M. The model is 
free of discrete gauge anomalies. In addition to giving rise to realistic masses and mixing 
patterns for the leptons and quarks, a unique feature of the SU(5) x T' model is that the CP 
violation, predicted by the model, is entirely geometrical in origin. This interesting aspect of 
the model considered is a consequence of the special properties of the group T' . The model 
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was shown to predict approximately tri-bimaximal neutrino (TBM) mixing with deviations 
which generate a non-zero sin 6*13 = sin 6 C / (3v2), where 6 C is the Cabibbo angle, and a Dirac 
phase 5 = 221° = 57r/4. A rather detailed description of the lepton sector of the model is 
given in Section 2. 

We have derived first the predictions of the model for the light Majorana neutrino masses, 
rrii, i = 1,2,3. We have found that in the model considered only light neutrino mass 
spectrum with normal ordering (or normal hierarchy) is possible. In the standardly used 
convention this implies m\ < m>2 < m^. The masses rrii depend on two parameters of the 
model which, as we have shown, can be fixed using the values of Am^, which drives the solar 
neutrino oscillations, and the ratio r = Am^/Am^, Am^ being the neutrino mass squared 
difference responsible for the dominant atmospheric and oscillations. Given the fact 
that Am^ and r have been determined experimentally with rather small uncertainties, the 
light neutrino masses are predicted to lie in relatively narrow intervals around the values 
mi = 1.14 x 10~ 3 eV, m 2 = 8.78 x 1(T 3 eV, m 3 = 4.89 x 1(T 2 eV (Fig. [[]). As a consequence, 
the model provides specific predictions for the sum of the light neutrino masses in terms 
of Amgj and r: m 1 + m 2 + m 3 = Am| 1 (l + (20r - l) 2 + (20r - 3) 2 )(1 - (20r - l) 4 )^. 
Numerically, we have found: mi + m 2 + m 3 = 5.9 x 10 -2 eV. 

The model provides also specific prediction for the effective Majorana mass in neutrinoless 
double beta (( / 9 / 9) 0! y-decay), |<m>| . The two Majorana phases 0:21 and a^i in the neutrino 
mixing matrix (see eqs. (Q]- (El)) are predicted (to leading order) to have CP conserving 
values: (X21 = and 031 = tt. We found that the effective Majorana mass predicted by 
the model is given approximately by: |<m>| = \(2mi + 7712) /3 — v / 2(^2 — mi)(sin0 c )/9\. 
Numerically, we obtained |<m>| = 3.4 x 10~ 3 eV. 

In the last part of this work we have derived detailed predictions for the rates of the 
lepton flavour violating (LFV) charged lepton radiative decays \x — > e + j, t — > e + 'j and 
r — > /i+7. This was done within the commonly employed mSUGRA SUSY breaking scenario 
J43|. The values of the mSUGRA parameters tan/3 and Aq are chosen from the intervals 
tan f3 = 3-i- 50, A$ = -j- 7mp, which are compatible with the constraints obtained by the 
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5l| and CMS |52| experiments at LHC. The values of the other two relevant 



mSUGRA parameters m „ and ,n m , are chosen from interval favored by the g.oba. data 
analysis performed in [53|. The data set used in this analysis includes in addition to the 



results of the ATLAS and the CMS experiments, the data on the muon (g — 2), on the 
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precision electroweak observables, on #-physics observables, astrophysical data on the cold 
dark matter density, as well as the limits from the direct searches for Higgs boson and 
sparticles at LEP. Based on the results obtained in 53], m and m^j are chosen from, or 



to vary in, the intervals 50 GeV < m < 400 GeV and 300 GeV < m 1/2 < 800) GeV. The 
branching ratios are calculated for masses of the heavy Major ana neutrinos = M = 10 12 
GeV. The GUT scale used is Mx = 2 x 10 16 GeV. The results of the this analysis are reported 
graphically in Figs. H] - [TO] 

One specific prediction of the SU(5) x T' model of flavour considered is that the quantity 
|(y^Yj) Zi /| on which the I — > /' + 7 decay branching ratios depend, I — fi, r, I' — e, fi 
(mi > mi/), Y v being the matrix of neutrino Yukawa couplings in the basis in which the 
RH neutrino Majorana mass term and the matrix of charged lepton Yukawa couplings are 
diagonal, is a function of the PMNS neutrino mixing matrix U and neutrino masses rrii only, 
KY^Yj)^! = \UijmjUj*,j\. Thus, in the model considered and for given I and /' this quantity is 
fixed numerically with small uncertainties. Together with the fact that the heavy Majorana 
neutrinos are degenerate in mass, this feature of the model implies that the double ratios 
of the branching ratios #(21/31) = (##(/x — > e + j)/BR(t — > e + j))BR(t — > ez/ T z/ e ) and 
#(21/32) = (BR(fi -> e + j)/BR(t fi + j))BR(t ev T v e ) do not depend not only 
on the mSUGRA parameters, but also on the heavy Majorana neutrino masses and the 
GUT scale M x : i?(21/31) = \(Y u Y})^ e \ 2 /\{Y v Y}) re \\ #(21/32) \{Y v Y})^ e \ 2 /\(Y u Y}) Til \ 2 . 
Numerically we have found: #(21/31) = 0.21, #(21/32) ^ 7.4 x 10" 3 . These values are 
one of the characteristic predictions of the model with SU(5) x X" symmetry considered. 
In particular, the r — > fi + 7 decay branching ratio, BR(r — > \i + j)/BR(t — > eu T u e ), is 
predicted to be bigger than the — > e + 7 decay branching ratio, BR([i — > e + 7), by a 
factor of 1.36 x 10 2 . We have found also that in a relatively large part of the mSUGRA 
parameter space considered and for M = 10 12 GeV, the // — > cy decay branching ratio can 
satisfy the MEG upper bound BR(fi — > e + 7) < 2.4 x 10~ 12 , but still can have a value in 
the range of sensitivity of the MEG experiment, BR( u — > e + 7) ^ 10~ 13 . 

The model with SU(5) x X" symmetry proposed in 16j,|l7| and investigated in the present 
article provides a unified description of the masses and mixing of the quarks and leptons, 
neutrinos included, as well as of the CP violation both in the quark and lepton sectors. 
This makes it a viable model of flavour, which possesses a number of appealing features. In 
the lepton sector the model provides specific predictions for the values of the light neutrino 



27 



masses, the neutrino mass spectrum, the values of the neutrino mixing angles, including the 
smallest one 6*13, and the Dirac and Majorana CP violation phases in the neutrino mixing 
matrix. All these predictions can and will be tested in the currently operating and future 
neutrino experiments. We are looking forward to the outcome of these tests. 
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